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Abstract 

We study gravitating monopoles and non-abelian black holes of SU(2) Einstein- 
Yang-Mills-Higgs theory coupled to a massless dilaton. The domain of existence 
of these solutions is bounded and decreases with increasing dilaton coupling 
strength. The critical solutions of this system are Einstein-Maxwell-dilaton 
solutions. 
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1 Introduction 



In SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory, with the Higgs field in the ad- 
joint representation, globally regular gravitating magnetic monopole solutions 0, §, 
emerge smoothly from the flat space 't Hooft-Polyakov monopole 0] , when gravity is 
coupled. Since the mass of the non-abelian monopole in flat space is proportional to 
the Higgs vaccum expectation value v, while the size of the monopole core is inversely 
proportional to v, the monopole size should become comparable to the Schwarzschild 
radius, when v becomes sufficiently large, and the non-abelian monopole should im- 
plode to a black hole ||. Thus regular magnetic monopoles should exist only in a 
limited domain of the EYMH parameter space. 

Indeed, as revealed by numerical analysis [I], [2|. |3|] , regular gravitating monopoles 
exist only up to a maximal value a max of the coupling constant a (which is propor- 
tional the Higgs vacuum expectation value v and the inverse Planck mass). Beyond 
ctmax the only solutions with unit magnetic charge are embedded abelian solutions, 
namely Reissner-Nordstr0m (RN) black hole solutions. Thus a max limits the domain 
of the EYMH parameter space, where non-abelian monopole solutions exist. 

Besides the maximal value of a, a max , the critical value of of particular 

importance. As a approaches the critical value a CT < a max the non-abelian 
monopole solutions approach a critical solution with a degenerate horizon. Outside 
the horizon the critical solution corresponds to an extremal RN solution with unit 
charge, while it is non-singular inside the horizon. Thus at a CT , the non-abelian 
monopole branch bifurcates with the branch of extremal RN black holes 0, |3| . 

Distinct from embedded RN black holes, genuine non-abelian magnetically charged 
black hole solutions emerge from the globally regular magnetic monopole solutions, 
when a finite regular event horizon is imposed [p], |2j, [Tj. Representing counterexam- 
ples to the "no-hair" conjecture for black holes, these black hole solutions may be 
viewed as "black holes within monopoles" 0, and have been interpreted as bound 
monopole-black hole systems M. 

The domain of existence in parameter space of genuine non-abelian black hole 
solutions is also bounded [|], ||. In particular, for a fixed value of the coupling 
constant a (a < a max ), the domain of existence of non-abelian black hole solutions is 
bounded by a maximal value of the horizon radius x/^max (which depends on a), and 
only relatively small non-abelian black hole solutions exist, as predicted on general 
grounds ||. 

Here we investigate the influence of a massless dilaton on the gravitating non- 
abelian monopole solutions, since dilatons appear naturally in many unified theo- 
ries, including string theory. The influence of a massless dilaton on the non-abelian 
monopole solutions in flat space has remarkable similarities to the influence of gravity 
on the monopole solutions HT0| . In particular, in flat space the non-abelian dilatonic 
monopole solutions exist only below a maximal value of the coupling constant 7 
(which is proportional the Higgs vacuum expectation value v and the dilaton cou- 
pling constant). As 7 approaches a critical value, the non-abelian dilatonic monopole 
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solutions approach a critical solution, which corresponds to an embedded abelian 
dilatonic monopole solution. Thus the branch of regular non-abelian monopole solu- 
tutions bifurcates with the branch of singular abelian monopole solutions . 

Considering the influence of both gravity and a dilaton, we expect that the do- 
main of existence in parameter space of the non-abelian monopoles in Einstein- Yang- 
Mills-Higgs-dilaton (EYMHD) theory should also be bounded. This expectation is in 
distinction to arguments from string theory, which suggest, that regular non-abelian 



dilatonic monopole solutions exist for any value of the coupling constant [11]. We fur- 



thermore expect that the non-abelian dilatonic monopole solutions should also show 
a critical behaviour, when the respective coupling constants are varied. In particu- 
lar, the branch of gravitating non-abelian dilatonic monopoles should bifurcate with 
the corresponding branch of embedded abelian solutions, namely Einstein-Maxwell- 



dilaton (EMD) solutions p 



We also investigate the influence of a dilaton on the non-abelian black hole solu- 
tions, and we determine their domain of existence in parameter space. Furthermore, 
we obtain simple relations between the dilaton field and the metric, and the dilaton 



charge and the mass for the monopoles and black holes [13]. 

We review the EYMHD action in section II, where we discuss the spherically 
symmetric ansatz, the field equations, the boundary conditions and the relations 
between metric and dilaton field. We discuss our numerical results for the globally 
regular monopole solutions in section III and for the non-abelian black hole solutions 
in section IV. The conclusions are presented in section V. We briefly review the 
numerical procedure in Appendix A. 



2 SU(2) Einstein- Yang- Mills- Higgs- Dilaton Equa- 
tions of Motion 

2.1 SU(2) Einstein- Yang-Mills-Higgs-dilaton action 

We consider the action of the EYMHD theory 

S = S G + S M = J L Gv ^gd 4 x + J L M ^gd 4 x . (1) 

The gravity Lagrangian Lq is given by 

^L G = , (2) 

where G is Newton's constant. The matter Lagrangian is given in terms of the gauge 
fields A/, the Higgs fields $ a , and the dilaton field * (a = 1, 2, 3) 

4ttL a/ = --e 2 ** F£,F' H '> a - -d^W - -D^ a D^ a - e - 2K ^V(^ a ) , (3) 
1 2 2 
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with Higgs potential 

V($ a ) = j(<f> a $ a -v 2 ) 2 , (4) 
The non-abelian field strength tensor is given by 

F; v = dpAZ - d u A; + ge ahc A\A c u , (5) 

and the covariant derivative of the Higgs field in the adjoint representation by 

D^ a = d^ a + ge abc A b ^ c , (6) 

g denotes the gauge field coupling constant, k the dilaton coupling constant, A the 
Higgs field coupling constant, and v the vacuum expectation value of the Higgs field. 

2.2 Spherically symmetric ansatz 

To construct static spherically symmetric dilatonic monopoles and black holes, we 
employ for the metric the static spherically symmetric Ansatz with Schwarzschild- 
like coordinates 

ds 2 = -A 2 (r)N(r)dt 2 + N-\r)dr 2 + r 2 d9 2 + r 2 sin 2 9d 2 V , (7) 
with the metric functions A(r) and 

N(r) = 1 - , (8 ) 

where the metric function m{r) is the mass function. For the gauge and Higgs fields, 
we use the purely magnetic hedgehog ansatz j|, [l], || 

Af = A t a = , (9) 

Af = - l ~ K ^ e v a , V= 1- ^ (r) singeg" , (10) 

$ a = v H(r)e r a . (11) 
A non-vanishing time-component of the gauge field would lead to dyon solutions 



14, 15 



The dilaton is a scalar field depending only on r 

* = #(r) . (12) 
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2.3 Field equations 

Let us now introduce dimensionless coordinates and fields 

$ -qj 

x = rgv , a = mgv , <p = — , ib — — . (13) 

v v 

The Lagrangian and the resulting set of differential equations then depend only on 
three dimensionless couplings constants, a, [3 and 7, 

r- M w V\ M H kM w 

a = VGv = — — , = = , 7 = kv = , 14 

gMpi g y/2M w g 

where M\y = gv, M# = \/2Xv and M-p\ = 1/y/G. Special cases are k = y/G 
(i.e. 7 = a), corresponding to string theory, and n 2 = Qi2±pl ([ e ^ — ^2±n a ^ 
corresponding to (4 + n)-dimensional Kaluza-Klein theory |T2| . 

Inserting the Ansatz into the Lagrangian and varying with respect to the matter 
fields yields the Euler-Lagrange equations 

(e 2 ^ANK'Y = A(e 2 ^ K<yR2 ~ ^ + H 2 K) , (15) 

x 2 

(x 2 ANH')' = AH(2K 2 + f3 2 x 2 e~ 2 ^(H 2 - 1)) , (16) 

(x 2 AN^)' = 2 1 A(e 2 ^(N(K') 2 + ^ '"^ - ) - e ~ 2 ^^(H 2 - l) 2 ) , (17) 

where the prime denotes the derivative with respect to x. 
Variation with respect to the metric yields the Einstein equations 

Gpv = 2a% u , (18) 

where T^ u is the energy-momentum tensor 

4vrT^ = g, u L M - 2^ . (19) 

With the following combinations of the Einstein equations 

Gtt = 2a 2 T tt = -2a 2 A 2 NL M , (20) 

g"G rr - g*G u = "4« 2 A^ , (21) 
we obtain for the metric functions /i and A the differential equations 

J = a 2 ( e 2 ^N(K') 2 + -Nx 2 (H') 2 + -±-{K 2 - ife 2 ^ + K 2 H 2 
V 2 2ar 

+ ^x 2 (H 2 - 1)V 2 ^ + \nx 2 ^) 2 ^ , (22) 

A' = a 2 xA (^—^p— e2 ^ + (H') 2 + (4>') 2 ^J ■ (23) 

For 7 = 0, the dilaton decouples, ip = 0, and the set of equations reduces to the 
EYMH equations studied previously [IT], 0, 0. 



2.4 Boundary conditions 

For the system of three second order equations, Eqs. (|15|) -( fL7]) , and two first order 
equations, Eqs. fl22|) -(E?3|), we need to specify eight boundary conditions. 



2.4.1 Dilatonic monopoles 

Let us first consider the boundary conditions for static globally regular solutions with 
finite energy, which are asymptotically flat. Regularity at the origin requires the 
boundary conditions 

K(0) = 1, H(0) = 0, a«VU=o = 0, M0) = 0. (24) 

At infinity asymptotic flatness and finite energy impose the boundary conditions 

K{oo) = , H(oo) = 1 , V(oo) = , A(oo) = 1 , (25) 



where the condition on ip can be imposed because of dilatational symmetry ||16|| , and 
the condition on A fixes the time coordinate. Note, that the value of the metric 
function /i at infinity determines the dimensionless mass of the solutions, /i(oo)/a 2 . 



2.4.2 Dilatonic black holes 

To construct black hole solutions which possess a regular event horizon at radius Xh, 
corresponding boundary conditions must be imposed at the horizon. Regularity at 
the horizon requires for the matter functions 

N'K' = K{R2 ~ 1} + e~ 2 ^H 2 K , (26) 

ar 

N'H' = ^(2K 2 + (3 2 x 2 e~ 2 ^(H 2 - 1)) , (27) 

jvy = ^(e 2 ^(N(K') 2 + ( ^ 2 ~ 1)2 ) - e~ 2 ^^(H 2 - l) 2 ) , (28) 



x 



2x 2 ' 4 



and for the metric functions 

MK) = y , (29) 

and A(xh) < oo. At infinity, the black hole solutions satisfy the same boundary 
conditions as the globally regular solutions, Eqs. (|25|) . 



2.5 Relations between metric functions and dilaton field 

In Einstein- Yang-Mills-dilaton (EYMD) theory, relations between the metric func- 
tions and dilaton field are known |13|]. In the following we demonstrate, that corre- 
sponding relations hold in EYMHD theory. 
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Let us introduce the abbreviations 

"2 _ 

~2x 2 ' I 



E A = N{K'f+ {K \ 2 1)2 , Fy = f3 2 -r(H 2 — I) 2 , (30) 



F# = ^iVa; 2 (iF) 2 + K 2 H 2 , F = e 2 ^ . (31) 
In terms of these the dilaton equation and the first Einstein equation read, 

(x 2 ANip')' = 2 1 A(FE A - F^Ey) , (32) 

H' = a 2 {FE A + E H + F~ l E v + ijVx 2 (V') 2 ) . (33) 
Contraction of the Einstein equations yields for the curvature scalar R 

- R = 2a 2 (-iV(V') 2 - \f~ 1 E v - ^E H ) . (34) 

Combination of these equations gives |T3| 

(x 2 AN^)' = - \x 2 R) . (35) 

With 

Ax 2 R = - [x 2 (2A'N + AN')] ' + 4A// , (36) 
we obtain for the dilaton field the equation 

(ANx 2 ip'y = (x 2 {2A'N + AN'))' . (37) 

This equation can be integrated, yielding 

^ (ln(A 2 A0V + ^ 

2a 2 \ V ') ANx 2 

with integration constant C. 

We can now derive the global relations between mass /j,(oo)/a 2 and dilaton charge 
D, defined via 

M x ) *_ztf? . (39) 
x 

Considering first the regular solutions, we integrate Eq. ( P7[ ) from zero to infinity, and 
take into account the asymptotic behaviour of the dilaton field, Eq. (|3"5|), and of the 
metric functions 

A'(x) oc 0(x~ 3 ), N'(x) oc + 0(x~ 3 ) , (40) 

x l 

as well as N(oo) = A(oo) = 1. This yields the relation between the dilaton charge 



and the mass of the solution 13 



D = ^. (41) 

cr 



Consequently the integration constant C in Eq. ( PB|) vanishes for regular solutions. 
Integrating Eq. (|38|) from x to infinity then gives a simple relation between dilaton 
field and tt-component of the metric [13| 



rir) =-L-ln(A 2 N) = ^\n(-g tt ) . (42) 



For black hole solutions we integrate Eq. (p7|) from the horizon to infinity. With 
the appropriate boundary conditions we thus obtain for black holes the relation |13| 

7 1 2 



D = ^ (M°o) - -xWLj . (43) 
The integration constant C therefore does not vanish for black holes. 



3 Dilatonic Monopoles 

The non-abelian magnetic monopole solutions of Yang-Mills-Higgs (YMH) theory [|J 
persist in EYMH theory [|l], |2|, [3], when gravity is coupled, and in YMHD theory [ [TCfl , 
when a dilaton field is coupled, for small values of the respective coupling constants. 
In both cases, the solutions bifurcate with an abelian solution at a critical value of 
the respective coupling constant. 

Here we construct non-abelian magnetic monopoles, coupled to both gravity and 
a dilaton field. We solve the set of differential equations Eqs. (|15D-(|17D and Eqs. (^)- 
(|23|) numerically (see Appendix A for details), subject to the boundary conditions 
Eqs. (^4[)-(^5l). We show, that the gravitating dilatonic non-abelian magnetic monopole 
solutions also bifurcate with an abelian solution, an extremal EMD solution |12] . 



3.1 Monopoles in EYMH theory 

Let us begin by briefly recalling the EYMH results for gravitating monopoles obtained 
in [|1], 0, [3], corresponding to the limit 7 = of EYMHD theory. 

For a fixed value of the Higgs self-coupling (3, a branch of gravitating non-abelian 
monopole solutions emerges smoothly from the flat space monopole solution, when 
a is increased. This branch exists on the interval < a < a max . We refer to this 
branch as main branch. The maximal value of a, a max , decreases with increasing (3. 

For vanishing and small values of /3, a short second branch of non-abelian monopole 
solutions is present on the interval a cr < a < a max , and bifurcates with the main 
branch at a max . We refer to this branch as secondary branch. At a cr the secondary 
branch bifurcates with a branch of embedded abelian solutions, the extremal RN 
solutions of Einstein-Maxwell (EM) theory, carrying unit magnetic charge and mass 
/iRN(oo)/a 2 = 1/a. For larger values of (3, the secondary branch disappears, and 
« cr = «max 0! 01 ■ Increasing (3 further leads to a new phenomenon || [T7| . 

Let us now consider, how the solutions approach the critical solution as a — > a CT , 
while (3 is kept fixed. The metric function N(x) possesses a minimum, which decreases 
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monotonically as a — > a cr . At a CT the minimum of N(x) reaches zero at the value 
of the radial coordinate x = Xh = a cr . For x > Xh, the metric function N(x) of 
the critical solution is identical to the metric function N^(x) = (1 — a CT /x) 2 of the 
extremal RN solution, whereas it is different from the abelian function for x < Xh- 
Xh represents the event horizon of the RN black hole. 

Likewise, the matter field functions K(x) and H(x) approach the constant values 
Krn(x) = and H RN (x) = 1 of the extremal RN solution as a — > a CT for x > Xh, 
whereas they tend to non-trivial functions for x < Xh- Thus the critical solution 
reached at a cr is identical to an extremal RN solution outside the radius Xh, while it 
retains distinct non-abelian features inside the radius Xh- 



3.2 Monopoles in YMHD theory 

Let us now recall the YMHD results for dilatonic non-abelian monopoles in flat space 
flQfl , corresponding to the limit a = of EYMHD theory. The dilatonic non-abelian 
monopoles in YMHD theory have many features in common with the non-abelian 
monopoles in curved space. 

When the dilaton is coupled, the main branch of dilatonic non-abelian monopoles 
emerges smoothly from the 't Hooft-Polyakov monopole, when 7 is increased, while 
the Higgs self-coupling (3 is kept fixed. The main branch extends up to a maximal 
value of 7, 7 m ax- For vanishing and small values of (3, a secondary branch exists on 
the interval 7 cr < 7 < 7 max , and bifurcates with the main branch at 7 max . At ■y CT the 
secondary branch bifurcates with a branch of embedded abelian monopole solutions. 
For larger values of /3, the secondary branch disappears, and 7 cr = 7 max flQfl . 

Intriguingly, even the numerical value for the maximal coupling constant 7 max is 
almost identical to the maximal coupling constant ct max (for (3 = 0). 



3.3 Monopoles in EYMHD theory 

Let us now consider how the presence of both gravity and a massless dilaton field 
affects the non-abelian monopoles. Generalizing the observations from EYMH and 
YMHD theory, we expect, that non-abelian monopoles in EYMHD theory also show 
a critical behaviour, when the respective coupling constants a and/or 7 are varied. 
In particular, the emerging branch of gravitating non-abelian dilatonic monopoles 
should bifurcate with the corresponding branch of embedded abelian solutions. Here 
these should be Einstein-Maxwell-dilaton (EMD) solutions fll2| . Thus the domain 



of existence in parameter space of the gravitating dilatonic non-abelian monopoles 
should be bounded. 



3.3.1 (3 = 

We first discuss the non-abelian monopoles of EYMHD theory for (3 = 0. 

Starting from the flat space dilatonic monopole solution for some value of the 
coupling constant 7, 7 < 7 cr [jlOfl , we couple gravity by choosing a small but finite 
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value of the coupling constant a. Increasing a, we obtain the main branch of dilatonic 
non-abelian monopole solutions, which exists on the interval < a < a max . The short 
secondary branch of dilatonic non-abelian monopole solutions, present on the interval 
«cr < ot < ctmax, bifurcates with the main branch at a max . At a CT the secondary branch 
bifurcates with a branch of embedded abelian solutions, which correspond to extremal 
EMD solutions, carrying unit magnetic charge and mass 

/^emd(oo) _ 1 (Ad , 



This is in complete accordance with our expectation |T8 . 

Let us now consider in more detail how the non-abelian EYMHD monopole solu- 
tions approach the critical value of a, when 7 is held fixed (and (3 = 0). As in the 
EYMH system, the metric function N(x) of the EYMHD system develops a minimum 
that decreases monotonically along the main and secondary branch. However, along 
the secondary branch the area close to the minimum starts to form a plateau, which 
broadens as a — > a cr . At the same time the function N{x) becomes increasingly steep 
close to the origin. In the limit a — > a cr , the plateau extends up to x = 0, and the 
non-abelian function N(x) converges to the function Nemd( x ) of the extremal EMD 



solution with unit magnetic charge [12 1, and with the same values of a and 7. But 



since the value of the extremal EMD solution at the origin is given by |T2 

AWd(O) = , (45) 

while the boundary conditions for the non-abelian solution require N(0) = 1, con- 
vergence is pointwise. In contrast to the EYMH case, the minimum of the function 
N(x) thus does not reach zero in the limit a — > a cr , when a dilaton field is present, 
so no horizon forms. 

We illustrate the dependence of the metric function N(x) of the EYMHD monopoles 
on the coupling strength a for 7 = 0.5 in Fig. la. The metric function N(x) is shown 
for the maximal value of a, a max = 1.30702, and for a = 1.29459, a value very close to 
the critical value a cr , where a pronounced plateau is already visible. For comparison, 
the extremal EMD solution for a = 1.29459 is also shown. 

The matter functions K(x) and H(x) also converge to their constant EMD values, 
Kemd(%) = and H EM d(x) = 1, in the limit a — > a cr , Because of their respective 
boundary conditions at the origin, K(0) = 1 and H(0) = 0, convergence to the EMD 
solution is again pointwise. 

Likewise, the dilaton function ip(x) approaches the EMD dilaton function 

7 X_ 

^emdO^) = -5- — J K 1 _ T > ( 46 ) 
cr + 7^ A 

with 

x/a = X(l - ^-f^ , X_ = ( ^ . (47) 



X ' \ a 
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Since X = X_ corresponds to the origin x — 0, 4>emd(%) tends logarithmically to 
minus infinity for x — > 0. We illustrate the dilaton function ip(x) for 7 = 0.5 in 
Fig. lb for a max = 1.30702 and a = 1.29459, i. e. for the same parameter values as 
in Fig. la. Convergence of the non-abelian functions ip(x) to the abelian function 
^emd(^), also shown in Fig. lb, is clearly seen. In particular, the value of the dilaton 
function at the origin, ip(0), diverges for a — > a cr , in accordance with Eqs. (^)-([T7|). 

For a fixed value of 7, the mass fi(oo)/a 2 of the non-abelian solutions decreases 
with increasing a. In Fig. 2 we exhibit the dependence of the mass of the non-abelian 
monopoles on the coupling constant a. In particular, we exhibit the mass of the 
EYMH solutions for 7 = 00 and the mass of the EYMHD solutions for 7 = 1.0 and 
1.2. Also shown is the mass of the respective abelian solutions, namely the extremal 
RN solutions for 7 = and the extremal EMD solutions for finite 7. The figure nicely 
illustrates the bifurcation of the non-abelian solutions with the corresponding abelian 
solutions. (The short secondary branches are very close to the abelian branches.) 

Fig. 2 also shows, that for a fixed value of a, the mass fi(oo)/a 2 of the non-abelian 
solutions decreases with increasing 7. Interestingly, the mass of the critical solution, 
where the bifurcation takes place appears to be rather constant, as seen in Fig. 2. 
Given by Eq. (|44|), the mass of the critical solution depends only on the combination 
a 2 + 7 2 . Since a max depends on 7, we conclude, that 

(cw(7)) 2 + 7 2 -const . (48) 
As demonstrated in Fig. 3, the values a max (7) lie indeed approximately on a circle, 



a max + 7 2 ~ 1-4, which forms the boundary of the domain of existence of the non- 
abelian EYMHD monopole solutions. 

When varying 7 while keeping a fixed, we obtain consistent results. In particular, 
for fixed a, the main branch of dilatonic non-abelian monopole solutions exists on 
the interval < 7 < 7 max , while the short secondary branch is present on the interval 
7cr < 7 < 7max- At 7 cr the secondary branch bifurcates with the corresponding branch 
of extremal EMD solutions Jl9| . 

Thus the dilaton field affects the gravitating non-abelian monopoles in a number 
of ways. Varying 7 along the main and secondary branch while keeping a fixed, 
the non-abelian core of the monopole solution shrinks and at the same time, the 
mass decreases. At the critical value of 7 the core shrinks to zero size, and the mass 
reaches the value of an extremal EMD solution, representing a naked singularity. This 
critical behaviour differs distinctly from the EYMH case, where the exterior part of 
the critical solution coincides with an abelian extremal black hole, and the interior 
part of the critical solution retains the non-abelian features. 

One should keep in mind the scale invariance of EYMHD theory, though. The 
above discussion of the critical behaviour corresponds to the choice ip(oc) = for the 
value of the dilaton field at infinity, Eq. (|2~5"D . This implies that the value of the dilaton 
field at the origin diverges for 7 — > 7 cr , ^((J) — > —00. A rescaling of the solutions 



16| such that ^(0) = should yield a critical solution with a regular interior and 



a finite size monopole core, but with an asymptotically diverging dilaton field and a 
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diverging mass |TI| . 



3.3.2 (3 / 

We have performed a detailed analysis of the EYMHD monopole solutions for small 
values of (3. Here the features observed for (3 = persist. The domain of existence 
shrinks with increasing (3. But its boundary remains close to a circle for small values 
of P, becoming slightly deformed with increasing (3. This is illustrated in Fig. 3 for 

/3 = 1 and (3 = 2. 

While we have made an extensive analysis of the EYMHD monopole solutions 
for vanishing and small (3, we have not been able to cover much of the parameter 
space for larger values of (3. When trying to determine the domain of existence in the 



a-7-plane for larger values of /3, we have encountered numerical problems f21| . Here 



new phenomena appear [p2]1 , which make the numerical analysis extremely difficult. 



Certainly further investigation of the non-abelian EYMHD monopole solutions for 
large values of j3 is called for, but employment of another numerical procedure might 
be necessary. 



4 Dilatonic Black Holes 

Non-abelian black holes in EYMHD theory have many features in common with non- 
abelian black holes in EYMH theory. In particular, they are limited in size and they 
tend to critical abelian solutions at critical values of the horizon radius. 

The black hole solutions of EYMHD theory depend on four parameters, a, [3, 7 
and the horizon radius x\ t . Here we limit the discussion to (3 = 0. We construct 
the black hole solutions, by solving the set of differential equations Eqs. (|15l)-(|17l) 
and Eqs. (|22"1) -([23"1) numerically (see Appendix A for details), subject to the boundary 
conditions Eqs. (p5|)-(p^), thus the black hole solutions are asymptotically flat and 
they possess a regular horizon. 

4.1 Black holes in EYMH theory 

Let us again briefly recall some properties of the non-abelian black holes in EYMH 
theory fl], ||, 0, f| . For a fixed value of the coupling constant a (a < ct max ) , black hole 
solutions emerge from the globally regular solution, when a finite horizon radius is 
imposed. They persist up to a maximal value of the horizon radius £/i,max, depending 
on a. Thus the domain of existence of genuine non-abelian black hole solutions is 
bounded, and only relatively small non-abelian black hole solutions exist |J. 

We dinstiguish two regions of black holes solutions with different critical be- 
haviour. For a fixed value of a > Vo/2, the critical solution approached corresponds 
to an extremal RN solution. In contrast, for a fixed value of a < v3/2, the critical 
solution approached corresponds to a non-extremal RN solution. In particular, for 
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values of a not too close to a = y/3/2, two branches of solutions appear, a main 
branch < Xh < a^,max an d a secondary branch x^ ct < Xh < ^/i,max- 

4.2 Black holes in EYMHD theory 

Let us first consider the domain of existence of the non-abelian black hole solutions. 
We recall, that the domain of existence of the non-abelian EYMHD monopole solu- 
tions is approximately bounded by a circle in the a-7-plane. Imposing a small but 
fixed horizon radius Xh, the domain of existence of the non-abelian EYMHD black 
hole solutions decreases with respect to the domain of existence of the regular solu- 
tions. With increasing size of the horizon, the domain of existence of the black hole 
solutions decreases further. Since the maximal value of 7, 7 max (a = 0) decreases 
stronger with increasing Xh than the maximal value of a, a max (7 — 0), the domain of 
existence of the non-abelian black hole solutions is no longer approximately bounded 
by a circle in the a-7-plane for finite values of Xh- This is demonstrated in Fig. 4, 
which shows the domain of existence in the a-7-plane for the non-abelian black hole 
solutions with horizon radii Xh = 0.1 and Xh = 0.3, in addition to the domain of 
existence of the non-abelian monopole solutions. 

In the a-x^-plane the domain of existence of the non-abelian black hole solutions 
decreases with increasing 7, as demonstrated in Fig. 5 for 7 = 0, 7 = 0. 5 and 7 = 1.0 
PSfl . In the EYMH system, the domain of existence is separated into two distinct 
regions by the straight line a = Xh, representing the extremal RN solutions. In 
particular, the maximum of the boundary curve X/ ljmax (a) occurs on this line for 
a = y/3/2. As seen in Fig. 5, the boundary of the domain of existence in the a- 
2^-plane changes drastically with increasing 7. The boundary curve becomes rather 
smooth and the maximum disappears. We note, that in contrast to the EYMH case, 
in the EYMHD case abelian solutions exist for arbitrarily small values of the horizon 
radius Xh, since extremal EMD solutions have x^ = 0. 

While Fig. 5 demonstrates the dependence of the maximal value of a, a max , on 
Xh for fixed 7, we now consider the dependence of the critical value of a, a CT , on Xh 
for fixed 7. In the following we restrict the discussion to 7 = 0.5. 

For fixed Xh and increasing a, the main branch of solutions extends from a = 
to a max (xft). At a max (xh) the main branch bifurcates with the secondary branch. 
Depending on the value of Xh, we observe two qualitatively different types of behaviour 
for the secondary branch. When < 2/1(7 = 0-5) < 0.22, the secondary branch ends 
at a finite critical value of a, whereas when 0.22 < 2^(7 = 0.5) < 0.48, the secondary 
branch extends all the way to a = 0. This is demonstrated in Fig. 6, where the values 
of the matter functions at the horizon, K(xh) and H(xh), are shown for Xh = 0.2 and 
x h = 0.3. 

When decreasing a along the secondary branch for < 2^(7 = 0.5) < 0.22, it 
first appears, that the solution tends to an extremal EMD solution on the interval 
2 > Xh with a divergent derivative of the metric function iV at Xh- Indeed, the 
mass ratio /i(oo)//i EMD (oo), remains close to one on a large part of the secondary 
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branch, as seen in Fig. 7. Decreasing a further, however, indicates that the critical 
solution is not an extremal EMD solution, but a non-extremal EMD solution. While 
for a — > a cr the matter functions K and H approach their abelian values, K(xh) —> 
and H{xh) — > 1, as seen in Fig. 6 for Xh = 0.2, the mass ratio /i(oo)//i EMD (oo), 
tends to a value significantly greater than one, as seen in Fig. 7. The profiles of the 
metric function N and the dilaton function ip support the conclusion as well, that 
the critical solution is a non-extremal EMD solution. In particular, the plateau of 
the metric function N, that forms when a is decreased from a max , disappears again, 
when a is decreased further towards a cr . 

For 0.22 < Xhi'y = 0.5) < 0.48 no critical value of a exists, instead the secondary 
branch extends all the way to a = 0, where K(xh) and H(xh) reach intermediate 
values in the interval ]0 : 1[ in the limit a — > 0, as demonstrated in Fig. 6 for 
Xh = 0.3. The corresponding mass ratio A*(oo)///emd(oo) for Xh = 0.3 is shown in 
Fig. 7. 

5 Conclusions 

We have constructed non-abelian monopoles and black holes in EYMHD theory, for 
vanishing and small Higgs self-coupling. Like the non-abelian monopole solutions of 
EYMH theory, the non-abelian monopole solutions of EYMHD theory reach critical 
solutions, when a critical value of the Higgs vacuum expectation value is reached. 
Therefore, the domain of existence of non-abelian monopole solutions is bounded also 
in EYMHD theory. 

As the critical solution is approached, the branch of non-abelian EYMH monopole 
solutions bifurcates with a branch of embedded abelian solutions, corresponding to 
extremal RN solutions with unit charge. In the presence of a massless dilaton field, 
the branch of non-abelian EYMHD monopole solutions also bifurcates with a branch 
of embedded abelian solutions, which here correspond to extremal EMD solutions 
with unit charge. Thus without the dilaton field the critical solution is associated 
with a black hole solution, whereas in the presence of the dilaton field the critical 
solution is associated with a naked singularity. 

Like the non-abelian black hole solutions of EYMH theory, the non-abelian black 
hole solutions of EYMHD theory are limited in size. Only relatively small black hole 
solutions exist ||. The critical solutions reached at a critical value of the horizon 
radius Xh also are associated with embedded abelian solutions, corresponding to RN 
black hole solutions in EYMH theory and to EMD black hole solutions in the presence 
of a dilaton field. 

The boundary of the domain of existence of the non-abelian monopole solutions 
is approximately a circle in the a-7-plane. For non-abelian black hole solutions the 
domain of existence decreases with increasing horizon radius Xh in the a-7-plane. In 
the a-x/i-plane the domain of existence decreases with increasing 7, exhibiting a much 
smoother boundary line for finite values of 7 than for 7 = 0. 

Besides presenting numerical results we have also derived an analytical formula 
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relating the dilaton function and the metric in the case of the dilatonic monopoles, 
and we have given analytical relations between the mass and the dilaton charge for 
monopoles and black holes. 

We have not yet considered the radial excitations of the dilatonic monopoles and 
black holes, which should be present, in analogy to the radial excitations of the EYMH 
solutions |jj 0. 

EYMH theory also possesses static axially symmetric monopole and black hole 
solutions, which are not spherically symmetric These black hole solutions, in 

particular, show, that Israels theorem cannot be extended to EYMH theory [ ^5] . 
Currently, we are investigating such static axially symmetric solutions also in EYMHD 
theory. 

Acknowledgements: One of us (B.H.) acknowledges the Belgium F.N.R.S. for fi- 
nancial support. 



6 Appendix A: Numerical Procedure 

To construct the dilatonic monopole and black hole solutions numerically, we employ 
a collocation method for boundary-value ordinary differential equations developed 
by Ascher, Christiansen and Russell p6f. The set of non-linear coupled differential 



equations Eqs. ([TBI) - fll7D and Eqs. ([22)) - (|23|) is solved using the damped Newton 
method of quasi-linearization. At each iteration step a linearized problem is solved 
by using a spline collocation at Gaussian points. Since the Newton method works 
very well, when the initial approximate solution is close to the true solution, the 
dilatonic monopole and black hole solutions for varying parameters a, /3, 7 and Xh 
are obtained by continuation. 

The linearized problem is solved on a sequence of meshes until the required accu- 
racy is reached. For a particular mesh Xi = x\ < X2 < ... < xn+i = x a , where Xi and 
x are the boundaries of the interval, and hi = Xi+i—Xi, h = maxi<j<jv/ij, a collocation 
solution v (x) = (vi,v 2 , ...Vd) is determined. Each component v n (x) G C™"" 1 ^,,!^ 
is a polynomial of degree smaller than k + m n , where m n is the order of the n-th 
equation, and k is an integer bigger than the highest order of any of the differen- 
tial equations. The collocation solution is required to satisfy the set of differential 
equations at the k Gauss-Legendre points in each subinterval as well as the set of 
boundary conditions. 

When approximating the true solution u n (x) by the collocation solution v n (x), an 
error estimate in each subinterval x G [xi,Xi+i) is obtained from the expression 

_ „W(rM\,., - n ,\o,( k + m ")(^A\h'y +Tn ™- 1 + 0(h k+mn ~ l+1 ) (49) 



K'K X ) -K J K x )m = °n,l\K "MX 



I = 0, m n — 1 , n = l,...,d, (50) 

where c n ^ are known constants. Also, using Eq. (|50|) a redistribution of the mesh 
points is performed to roughly equidistribute the error. With this adaptive mesh 
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selection procedure, the equations are solved on a sequence of meshes until the suc- 
cessful stopping criterion is reached, where the deviation of the collocation solution 
from the true solution is below a prescribed error tolerance [^(J . 

For the numerical solutions, we typically specified the error tolerance in the range 
10 -4 — 10 -6 . The number of mesh points used in these calculations was typically about 
350. We calculated the dilatonic monopole solutions on the finite interval [0, x Q ] and 
the black hole solutions on the finite interval [xh, x ], where Xh represents the horizon 
of the black hole, with x a ranging from 10 2 to 10 4 . 
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Figure la: The metric function N(x) of the non-abelian monopole solutions is shown 
for 7 = 0.5 and two values of a, the maximal value of a, a max = 1.30702, and 
a = 1.29459, a value very close to the critical value a CT . For comparison, the metric 
function N EMD (x) of the extremal EMD solution for 7 = 0.5 and a = 1.29459 is also 
shown. 
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Figure lb: Same as Fig. la for the dilaton function ip(x). 
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Figure 2: The dependence of the mass fi(oo)/a 2 on a is shown for the non-abelian 
monopole solutions for three different values of 7, 7 = 0, 7 = 1 and 7 = 1.2. For 
comparison, the dependence of the mass on a is also shown for the corresponding 
abelian solutions (RN and EMD, respectively). 
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Figure 3: The domain of existence of the non-abelian monopole solutions in the a-7- 
plane is shown for three values of the Higgs self-coupling, (5 — 0, (5 — 1 and (5 — 2. 



22 




Figure 4: The domain of existence of the non-abelian black hole solutions in the a- 
7-plane is shown for x^ — 0.1 and Xh = 0.3. Also shown is the domain of existence of 
the non-abelian monopole solutions. 
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Figure 5: The domain of existence of the non-abelian black hole solutions in the a- 
rr^-plane is shown for three values of the dilaton coupling constant, 7 = 0, 7 = 0.5 
and 7=1. 
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Figure 6: The values of the matter functions at the horizon, K(xh) and H(xh), of the 
non-abelian black hole solutions are shown for 7 = 0.5 and for Xh = 0.2 and x^, = 0.3. 
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Figure 7: The mass ratio /i(oo)//iEMD(oo) of the non-abelian black hole solutions 
and the extremal EMD solutions is shown for 7 = 0.5 and for Xh = 0.2 and x^ = 0.3. 
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